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Abstract 

The eletromagnetic field in a hnear absorptive dielectric medium, 
is quantized in the framework of the damped polarization model. A 
Hamiltonian containing a reservoir with continuous degrees of free- 
dom, is proposed. The reservoir minimally interacts with the dielectric 
polarization and the electromagnetic field. The Lagevin-Schrodinger 
equation is obtained as the equation of motion of the polarization 
field. The radiation reaction electromagnetic field is considered. For 
a homogeneous medium, the equations of motion are solved using the 
Laplace transformation method. 



1 Introduction 



In a homogeneous and nondispersive medium, the photon is associated with 
only the transverse cmponents of the electromagnetic field, in contrast in an 
inhomogeneous nondispersive medium, the transverse and the longitudinal 
degrees of freedom are coupled. In this case the quantization of the electro- 
magnetic field can be accomplished by employing a generalized gauge [1,2], 
that is, V.{s{r)A) = 0, where e{r) is the space dependent dielectric function. 
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Generalization of this quantization to an anisotropic nondispersive medium 
is straightforward by using of the gauge Y.'j=i — [3]- 

The quantization in a dispersive and absorptive dielectric, represents one of 
the most and interesting problems in quantum optics, because it gives a rig- 
orous test of our understanding of the interaction of hght with matter. The 
dissipative nature of a medium is an immediate consequence of its dissper- 
sive character and vice versa according to the Kramers- Kronig relations. This 
means that the validities of the electromagnetic quantization in nondissipa- 
tive but dispersive media is restricted to some range of frequencies for which 
the imaginary part of frequency dependent dielectric function is negligible, 
otherwise there will be inconsistencies. 

In the scheme of Lenac [4], for dispersive and nonabsorptive dielectric media 
starting with the fundamental equations of motion , the medium is described 
by a dielectric function s{f,u>), without any restriction on its spatial behav- 
ior. In this scheme it is assumed that there is no losses in the system, so the 
dielectric function is real for the whole space. The procedure is based on an 
expansion of the total field in terms of the coTiplcd cigcnmodes, orthogonality 
relations are derived and equal-time commutation relations are discussed. 
Huttner and Barnett [5] have presented a canonical quantization for the elec- 
tromagnetic field inside the dispersive and absorptive dieletrics based on a 
microscopic model in which the medium is represented by a collection of the 
interacting matter fields. The absorptive character of the medium is modeled 
through the interaction of the matter fields with a reservoir consisting of a 
continuum of the Klein-Gordon fields. In their model, eigen-operators for 
the coupled systems are calculated and the electromagnetic field has been 
expressed in terms of them, the dielectric function is derived and it is shown 
to satisfy the Kramers- Kronig relations. 

Matloob [6] has quantized the electromagnetic field in a linear isotropic 
medium by associating a damped harmonic oscillator with each mode of the 
radiation field. A canonical approach has been used to quantizing a damped 
quantum oscillator. The conjugate momentum is defined and a quantum 
mechanical Hamiltonian is introduced. 

Gruner and Welsh [7] have given a quantization scheme for the radiation 
field in dispersive and absorptive linear dielectric by starting from the phe- 
nomenological Maxwell equations, there the properties of the dielectric are 
described by a permitivity consistent with the Kramers- Kronig relations, an 
expansion of the field operators is performed based on the Green function 
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of the classical Maxwell equations which preserves the equal-time canonical 
commutation relations. In particular, in frequency intervals with approxi- 
mately vanishing absorption, the concept of quantization through mode ex- 
pansion for dispersive dielectrics is recognized. The theory further reveals 
that weak absorption gives rise to space-dependent mode operators which 
spatially evolve according to the quantum Langevin equations in the space 
domain. 

Suttorp and Wubs [8] in the framework of a damped polarization model, 
have quantized the electromagnetic field in an absorptive medium with spa- 
tial dependendence of its parameters. They have solved equations of motion 
of the dielectric polarization and the eletromagnetic field by means of the 
Laplace transformation for positive and negative times. The operators that 
diagonalize the Hamiltonian are found as linear combinations of canonical 
variables with coefficients depending on the electric susceptibility and the 
dielectric Green function. The time dependence of the electromagnetic field 
and dielectric polarization are determined. Present authors [9] have given a 
model in which a linear absorptive dielectric is modeled by two independent 
reservoirs which are called electrical and magnetic reservoirs. In this model, 
the electrical and magnetic polarization densities, for the reservoirs, are de- 
fined. The electrical and magnetic polarization densities interact with the 
displacement vector field and the magnetic field respectively. Both structural 
and Maxwell equations are obtained as the Heisenberg equations of motion. 
There are some another approaches for quantizing the electromagnetic field, 
see for example [10-19]. 

In this paper in the framework of the damped polarization model, we in- 
troduce a method in which the electromagnetic and dielectric polarization 
fields, interact with a reservoir through a minimal coupling term. We obtain 
the Langevin-Schro dinger equation as the equation of motion of the dielec- 
tric polarization and calculate the radiation reaction electromagnetic field in 
terms of the polarization field. Finally we solve the equations of motion in a 
homogeneous dielectric using the Laplace transformation. 

2 Quantum dynamics 

We take a model in which the polarization of the dielectric and the electro- 
magnetic field, as quantum fields, interact with a reservoir through a minimal 
coupling iterm. In this model we can obtain the Maxwell equations and the 
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Langevin-Schrodinger equation as the equations of motion of the polarization 
density of the dielectric. The vector potential for the electromagnetic field in 
the coulomb gauge, V.A = 0, can be expanded in terms of the plane waves 



A{f,t) = y f d?k, ^-4-, — 



^ K^,(t)e^'=->at (t)e-^'='1e(fc,A), (1) 



where a;^ = c\k\ and Eq is transitivity of vacuum, e{k, A) are polarization unit 
vectors 

e{k,X) ■ e{k,X') = 6x\', 

k-e{k,X) = 0, A =1,2, (2) 

Canonical momentum density of the electromagnetic field is defined by 



Mr, t) = -^e, J ^'W^T^^^^a W^*" - 4,(t)e-^^1e(^, A). (3) 



A=l ' 



The creation and annihilation operators at^, Uj^^, in any instant of time, 
satisfy the comutation relations 

[a^,,it),al^Xt)] = Sxx'5{k~k'), (4) 

and their time dependence is to be determined from the total Hamiltonian. 
Comutation relations (jU lead to 

[Mf,t),7rF^{f\t)] = ih6i^{r-r'), (5) 

where 5^j{f—r') = J (Pke^^'^'^~'^'\5ij — -j^), is the transverse delta func- 
tion. The Hamiltonian of the quantum electromagnetic field is 



2eo 2/io x=i 



The existence of electrical field inside a dielectric medium, cause electrical 
charges to move from their stable positions. Let a(f^ be the density of the 
displaced charges and Y{f, t) be their displacement, then the polarization 
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density induced in the medium is F = aY. We can take F as a quantum 
field and expand it in terms of the plane waves 

nr,t) = E / ^^fe(t)e^^^-> 4(t)e-^1u.(fc), (7) 

where Uu{k) — e{k^v) for v — 1,2 and uz{k) — k — ^. The sum over 

|fe| 

V — 1,2, gives the transverse component of Y and = 3, is the longitudinal 
component. We can write the canonical momentum density of the quantum 
field Y as 

Mr,t) = -ij:J c^^^^^^fe(i)e''^-"-4(t)e-^1^.(^). (8) 
Imposing the commutation relations 

K^.W><vW] = ^-"^(^-^"''), (9) 

on d^^ and rft^, leads to the following commutation relations for vector fields 

Y and Try 

[y,(f,t),7ry,(P,t)] = zM,,-5(f- ?). (10) 

Let us assume a spring force — p{r)uiQ{r)Y , exerted on an element of the 
medium with volume (i^r centered in r, where p(r), is the mass density of 
the displaced charges, then the Hamiltonian of the quantum field Y can be 
written as 




If the damping forces together with the electrical force are exerted on the 
elements of the dielectric, then in the damped polarization model, it will 
be assumed that the damping forces are caused by a quantum field, which 
is called the reservoir or the environment. In this model, the equation of 
motion of the quantum field Y and the reservoir together with the Maxwell 
equations, can be obtained from the Heisenberg equations using the total 
Hamiltonian 

J ^ 2p + ^Pi^^l^^)^ + ^ ^ ^ ' 

(12) 
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where (aY)", is the longitudinal component of the polarization density and 
is the Hamiltonian of the reservoir or environment defined as 

HR = j2l d'k j d^qhuj^,[blik, q, t)K{k, q, t) + h (13) 

with the commutation relations 

[K{k, q, t),bl{k', q', t)] = 6,,,6{k - k')5{q- q'). (14) 

The operator R{f,t), plays a crucial rule in the interaction between the 
reservoir with the electromagnetic field and the quantum field Y, by definition 

R=j2 fd'k f -^[f{u^^,f^K{k,q,t)e'^^Ur{u^,,f^bl{k,q^ 
v=i-' \/(27r)3 

(15) 

The function f{u^,r), is called a coupling function which is position depen- 
dent in an inhomogeneous medium. According to the total Hamiltonian ()12|) . 
the reservoir interacts with both the electromagnetic field and the quantum 
field y, while in previous models, the reservoir was interacting only with the 
field of the polarization density of the dielectric. 



2.1 The equation of motion of the vector field Y 

Because the damping forces are exerted on the elements of the dielectric, the 
quantum field F, describes a dissipative quantum system. Such a system 
can be described in terms of the Langevin equation [20] which has a broad 
and general application. This description can be formulated using a coupling 
between the system and a quantum mechanical heat-bath. The Heisenberg 
equation for the dissipative system, takes the form of Langevin-Schrodinger 
equation. The equation of motion of the position operator x(t), is 

mf + f dt'fiit - t')S{t') + V{x) = F/v(t), (16) 
Jo 

where the coupling with the heat-bath is described by two terms, a term that 
describes the absorption of energy by the heat-bath which is characterized 
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by the memory function and a fluctuating term, characterized by the 
operator valued noise force -F/v(t). Both terms are necessary for a consistent 
description of a quantum damped system. In the present model, with the 
proposed Hamiltonian (jl2|) . we can obtain the Langevin-Schrodinger equa- 
tion for the quantum field Y by combining the Heisenberg equations of the 
reservoir and the field Y. The equations of motion for the fields Y and vfy, 
can be obtained from the Heisenberg equations 

— = i\H VI = -aA~ R 
dt ~ h} ' ^ ~ p 

^ = ^[H,ny] = -pu;o^(f)f - ^(aF)ll, (17) 

so 

d'^Y dA a, - m dR dR , , 

p— + p^UrlY = -a- - -(a.y)" = _, (18) 

where E-^ — — ~ ^"JJ^^ ; ^ire transverse and longitudinal compo- 

nents of the electrical field respectively. Using ((HI), the Heisenberg equation 
for bu{k,q,t) can be obtained as 

K(k,q,t) = ^[H,b^{k,q,t)] = -iu^K{k,q,t) + 

with the formal solution 

K{k,q,t) = e-'^'i'Kik.l^) + 



Substituting b^{k,q,t), from in (fTHj) . gives 



at' 



JO 

(21) 
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where 



-f{t ~t',f) = — duiu^ \f{f,u)\'^ COS u{t-t'), (22) 

is the memory function and describes the absorption of energy of the medium 
by the reservoir. The field ^{f, t) 

af,t) = ij2fSkl ^!£^^^4/(f,c.,^)6,(^,g,0)e-^^E*+^«^-^^- 
J J J (27r)3 



v=\ ■ 



rir, u^)bl{k, q, 0)e+'-s*-n«.(9l. (23) 

is the noise field associated with the absorption and has a zero expectation 
value in the eigenstates of the reservoir. The equation ()21|) . is called the 
Langevin-Schrodinger equation for the damped quantum field Y. The term 

/o dt''y{t — t', r)Y{f, t'), usually generates a damping force for Y, for example 
if we take 

l/MP - (24) 

we obtain 

pY + pul{f)Y + l3Y{r,t) = aE^{f,t) + aE\\{f,t)+^{f,t), (25) 

which has a damping term proportional to the velocity and ^(r, t), is the 
noise field (|^ with the coupling function 



2.2 The radiation reaction 

In QED, a charged particle in quantum vacuum interacts with the vacuum 
field and its own field known as the radiation reaction. In classical electro- 
dynamics, there is only the radiation reaction field that acts on a charged 
particle in the vacuum. The vacuum and radiation reaction fields have a 
fluctuation-dissipation connection and both are required for the consistency 
of QED. For example the stability of the ground state, atomic transitions 
and lamb shift can only be explained by taking into account both fields. If 
self reaction was alone the atomic ground state would not be stable. When 
a quantum mechanical system interacts with the vacuum quantum field, the 
coupled Heisenberg equations for both the system and the quantum vacuum 
field give us the radiation reaction field, for example it can be shown that 
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the radiation reaction for a charged harmonic oscillator is |^ [21]. One can 
find the Heisenberg equation for the annihilation operators a^^^ 



2(27r)%o^fc 

(26) 



with the formal solution 



(27) 

Substituting a^^i^) from (jTfjl in E-^ = — we obtain 

= E^ + E^^, (28) 

where 



E^ = ^Jd'k^^^^Ja^,,{0)e--^'^^^^^^^ (29) 

is the transverse vacuum field and 
1 



1 d^r'air') 









X {Y{r',t') -[k.Y{r',t')]k} (30) 
is the transverse radiation reaction electrical field [21]. 

2.3 The Maxwell equations 

The Maxwell equations can be obtained as the Heisenberg equations of A 
and ttf 

at h eq 

-gf = -^[H,n,] = {aY) (31) 



after eliminating TTp, 



.„ = (aY)- — (32) 

By defining the displacement vector field ID as D = = (aY)^ + EqE^ 
and S = V X A, the equation (jH^ can be rewritten as 

V X S = //o^ (33) 

which is the familiar form of the Maxwell equation in a non-magnetic dielec- 
tric medium. In the absence of the external electrical charge density, we have 
D\\ = EoEW + iaff = or EW = -i^^A^ 



3 Solution of the Heisenberg equations 

The equation of motion defined in the previous sections, can be solved by 
using the Laplace transformation method. For any time dependent operator 
g{t), the forward and the backward ((^•^^''^(t)), Laplace transformation, are 
defined respectively as 

roo 

g^is) = / dtg{t)e-'\ 
Jo 

roo 

g\s) = dtg{-t)e-^K (34) 

Now taking the Laplace transformation of the equation ()21|) 

aY^{f, s) = aY/j{f, s) + eox{r, s)E^{f, s), 

aY\f, s) = aY^if, s) + eox{r, s)E\r, s), (35) 

where 

^(t**, 5) = ~-- , T(^5 '^) ~ STTt? ( dio — — - — ^ — ^ — 

(36) 
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The function x(r, s), is the Laplace transformation of the electrical suscep- 
tibility and 

aF/(f,s) = — x(r,s){e-^(f,s)+psf(f,0) + 7fy(f,0)-aA(f,0) 

a 

- i?(f,o) + 7(f,s)y(f,o)}, 

aY^{r,s) = —x{f,s){^{f,s) + psY{r,0)-7rY{r,0) + aA{r,0) 
a 

+ ^(f,0) + 7(f,s)y(f,0)}, (37) 

are the forward and backward Laplace transformations of the noise po- 
larization densities respectively. Another main equation can be obtained 
by taking the time derivative of the equation (jH^ and using the relation 

d'^E d^Y -, ^ 

-Aio^o-^ = ^o^-Q^ + V X V X (38) 

Equation (jHHj) . can be written in terms of the Laplace transformed compo- 
nents 

V X V X E^(f, s) + eQpos^e{r, s)E^{f, s) = ,P {r, s), 

V X V X E\f, s) + £o/ios^£(r, s)E\r, s) = J^(f, s), (39) 

where e{f,s) = 1 + x{f,s), is the Laplace transformation of the electrical 
permeability and 

//(f, s) = V X V X A{f, 0) - 

fioSTVpi^f, 0) - fios{aY)\f, 0) + fiosaY{f, 0) - fios'^aY^{f, 0), 
/(f, s) = -V X V X A{f, 0) - 

HosnF{r,0) - fios{aY)\f,0) + fiosaY{f,0) - pos'^aY^{f,0), 

(40) 

are called the source noise densities. In reference [8], the wave equations 
(jH^ have been solved using the Green function method. In the following for 
simplicity, we solve the Heisenberg equations ()18|) and ()31|) for an homoge- 
neous dielectric, that is when p,uJo,a,f are independent of the position r. 
In the case of an homogeneous and linear dielectric, in the absence of exter- 
nal charge density, the polarization charge density is zero and therefore the 
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vector field aY, has only the transverse components. In this case, using the 
Laplace transformation of (j31|) and (j35p 

V X V X {f, s) + fioeos'^e{s)A-^{f, s) = (XosaYl{r, s) + iJ.oeose{s)A(r, 0) 
-fioo:Y{f, 0) + fiQiTFif, 0), 

V X V X A^(f, s) + /io£oS^e(s)A^(r, s) = -HosaY^{f, s) + Hoeose{s)A{f, 0) 
+fioaY{f, 0) - fionpif, 0). (41) 

The equations ()41|) . can be solved easily using the Fourier transformation 



2 r - 
A{r, s) = J2 



A=l ' 
2 



« A=l 



a-^ J ^2(2^)3 ^1 fc2 + ^^^^52^(5) 



a 

+ C.C, (42) 

where the upper and lower signs (±), give A-f{f, s) and A''{f, s), respectively. 
Now A{'r,t) for t > is the inverse Laplace transformation of 74-^(r, s) and 
A{f, —t) for t > 0, is the inverse Laplace transformation of A^{f, s), that is 

A{f, t) = SOR[e''P{f, s)], t> 0, 

A{f, t) = SOR[e-''A^{f, s)], t<0, (43) 

where SOR[f{s)], means the sum of the residues of a complex function 
f. If the residues of the complex functions ■^T^'^fc ipuifj±is)x(s) 

P+^^e^lieis) ' ^^^^ respect to s have negative real parts, as usually it is the case 
for an absorptive dielectric, then we find the following behaviour in t 00 

Air, t) = d^k I -£L E 0)xH.,) ^..^^u:,.^ 

+ C.C. (44) 



12 



We can find the polarization field Y{f, t) by substituting A^{r^ s) and A^(r, s) 
from (j42j) in (j35|) and then taking the inverse Laplace transformation of 

Yf{f, s) and y^(r, s). 

4 Conclusion 

By taking a reservoir with continuous degrees of freedom and a suitable 
Hamiltonian, we could give a consistent quantization scheme for electro- 
magnetic field in an absorptive dielectric medium, in the frame work of the 
damped polarization model. In this model, the damped polarization field, 
was satisfying the quantum Langevin equation containing the electrical field 
as the source term. Using the Laplace transformation, the quantum Langevin 
equation solved and the susceptibility of the dielectric obtained in terms of 
the coupling function between the reservoir and the polarization field. Equa- 
tions containing a noise term for the vector potential and the electrical field 
obtained and solved which led to explicit forms of the vector potential and 
the electrical field. 
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